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Abstract 

Using magnetically trapped atomic hydrogen as an example, we investigate 
the prospects of achieving Bose-Einstein condensation in a dilute Bose gas. 
We show that, if the gas is quenched sufficiently far into the critical region of 
the phase transition, the typical time scale for the nucleation of the condensate 
density is short and of 0(h~/kBT c ) . As a result we find that thermalizing elastic 
collisions act as a bottleneck for the formation of the condensate. In the case 
of doubly-polarized atomic hydrogen these occur much more frequently than 
the inelastic collisions leading to decay and we are lead to the conclusion that 
Bose-Einstein condensation can indeed be achieved within the lifetime of the 
gas. 

PACS numbers: 67.65. +z, 32.80.Pj, 64.60.Qb 
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I. INTRODUCTION 



In the last few years it has been clearly demonstrated that not only charged ions but 
also neutral atoms can be conveniently trapped and cooled by means of electro-magnetic 
fields. Although the physics of the various ingenious scenarios developed to accomplish 
this is already interesting in itself JIJ, the opportunities offered by an atomic gas sample at 
very low temperatures are exciting too. Examples in this respect are the performance of 
high-precision spectroscopy, the search for a violation of CP invariance by measuring the 
electric dipole moment of atomic cesium || , the construction of an improved time standard 
based on an atomic fountain ||, and the achievement of Bose-Einstein condensation in a 
weakly-interacting gas. 

In particular the last objective is an important motivation for studying cold atomic gases 
and has been pursued most vigorously with atomic hydrogen [fH|. However, it was recently 
proposed that also the alkali-metal vapors cesium || and lithium |7j are suitable candidates 
for the achievement of Bose-Einstein condensation. We will nevertheless concentrate here 
on atomic hydrogen, because it still seems to be the most promising system for the ob- 
servation of the phase transition in the near future. Moreover, it has the advantage that 
the atomic interaction potential is known to a high degree of accuracy. As a result we can 
have confidence in the fact that the scattering length is positive, which is required for the 
condensation to take place in the gaseous phase ||, and small enough to rigorously justify 
the approximations made in the following for the typical temperatures (T ~ 10 fiK) and 
densities {n ~ 1 • 10 14 cm" 3 ) envisaged in the experiments. 

Due to the spin of the electron and the proton, the ls-hyperfine manifold of atomic 
hydrogen consists of four states which are in order of increasing energy denoted by \a), \b), 
|c), and \d), respectively. Only the |c) and \d) states can be trapped in a static magnetic 
trap, because in a magnetic field they have predominantly an electron spin-up component 
and are therefore low-field seeking ||. Furthermore, if we load a trap with atoms in these 
two hyperfine states, the |c) state is rapidly depopulated as a result of the much larger 
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probability for collisional relaxation to the high- field seeking \a) and \b) states which are 
expelled from the trap. In this manner the system polarizes spontaneously and we obtain 
a gas of |o?)-state atoms, known as doubly spin-polarized atomic hydrogen since both the 
electron as well as the proton spin are directed along the magnetic field. Unfortunately, 
such a doubly-polarized hydrogen gas still decays due to the dipole interaction between the 
magnetic moments of the atoms. Although the time scale 7"j ne ; for this decay is much longer 
than the time scale for the depopulation of the |c) state mentioned above, it nevertheless 



limits the lifetime of the gas sample to the order of seconds for the densities of interest |T0 



Having filled the trap with doubly-polarized atoms, we must subsequently lower the 
temperature of the gas to accomplish Bose-Einstein condensation. At present it is believed 
that the most convenient way to achieve this is by means of conventional [11],|12 or light- 
induced |13 evaporative cooling. In both cases the idea is to remove, by lowering the well- 
depth or by photon absorption in the perimeter, the most energetic particles from the trap 
and thus to create momentarily a highly nonequilibrium energy distribution that will evolve 
into a new equilibrium distribution at a lower temperature. According to the quantum 
Boltzmann equation describing this process, a typical time scale for the evolution is the 
average time between two elastic collisions r e ; = l/n(va), with (va) the thermal average 
of the relative velocity v of two colliding atoms times their elastic cross section a. Clearly, 
r e i must be small compared to r ine i to ensure that thermal equilibrium is achieved within 
the lifetime of the system. As a result, the minimum temperature that can be reached 
by evaporative cooling is about 1 [iK and indeed below the critical temperature of atomic 
hydrogen at a density of 1 • 10 14 cm" 3 . 

The previous discussion appears to indicate that a typical time scale for the formation 
of the condensate is given by r e i. However, this is not correct because simple phase-space 
arguments show that a kinetic equation cannot lead to a macroscopic occupation of the 
one-particle ground state: Considering a homogeneous system of N bosons in a volume V, 
we find from the Boltzmann equation that the production rate of the condensate fraction is 
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d N, 



=C rM(l + JV ff ) (1) 

dt N i n V K 0) ' K J 

where Nq is the number of particles in the zero-momentum state and C is a constant of 0(1). 

Hence, in the thermodynamic limit (N, V —>■ oo in such a way that their ratio n = N/V 

remains fixed) a nonzero production rate is only possible if a condensate already exists 



14j| and we are forced to conclude that Bose-Einstein condensation cannot be achieved by 



evaporative cooling of the gas. 

II. NUCLEATION 

In the above argument we have only considered the effect of two-body collisions. It is 
therefore legitimate to suggest that perhaps three or more body collisions are required for 
the formation of the condensate, even though they are very improbable in a dilute gas |15| . 



However, we can easily show that the same argument also applies to these processes: For a m- 
body collision that produces one particle with zero momentum we have 2m — 2 independent 
momentum summations, leading to a factor of V 2m ~ 2 . Moreover, the transition matrix 
element is proportional to V ■ V~ m due to the integration over the center-of-mass coordinate 
and the normalization of the initial and final state wave functions, respectively. In total the 
production rate for the condensate fraction is thus proportional to V 2m ~ 2 {y 1 " rn ) 2 V~ 1 {l-\-N^) 
or V _1 (l + Nq), which again vanishes in the thermodynamic limit if there is no nucleus of 
the condensed phase. As expected, the contributions from collisions that produce more than 
one zero-momentum particle have additional factors of V~ x and vanish even more rapidly if 
V -> oo. 

Clearly, we have arrived at a nucleation problem for the achievement of Bose-Einstein 
condensation which seriously endangers the success of future experiments. Fortunately, we 
suspect that the line of reasoning presented above is not completely rigorous because oth- 
erwise it implies that also liquid helium cannot become superfluid, in evident disagreement 
with our experience. Indeed, by using a kinetic equation to discuss the time evolution of the 
gas we have in effect neglected the buildup of coherence which is crucial for the formation of 



the condensate. Our previous argument therefore only shows that by means of evaporative 
cooling the gas is quenched into the critical region on a time scale T e i, not that Bose-Einstein 
condensation is impossible. To discuss that point we need a different approach that accu- 
rately describes the time evolution of the system after the quenching by taking the buildup 
of coherence into account exactly. Such a nonequilibrium approach was recently developed 
on the basis of the Keldysh formalism and can, in the case of a dilute Bose gas, be seen as 
a generalization of the Landau theory of second-order phase transitions [16||. As a conse- 
quence it is useful to consider the Landau theory first. This leads to a better understanding 
of the more complicated nonequilibrium theory and ultimately of the physics involved in the 
nucleation of Bose-Einstein condensation. 



A. Landau theory 

As an introduction to the Landau theory of second-order phase transitions we use the 
example of a ferromagnetic material |I7||. To be more specific we consider a cubic lattice 
with spins Si at the sites {i}. The Hamiltonian is taken to be 

H = -JY,S l -S 3 , (2) 

(ij) 

where J is the exchange energy and the sum is only over nearest neighbors. For further 
convenience we also introduce the magnetization 

tf=4ss- (3) 

% 

Physically it is clear that this model has a phase transition at a critical temperature T c of 
0(J I fee). Above the critical temperature the thermal fluctuations randomize the direction of 
the spins and the system is in a disordered (paramagnetic) state having a vanishing average 
magnetization (M) . However, below the critical temperature the thermal fluctuations are 
not large enough to overcome the directional effect of the Hamiltonian and the spins favor an 
ordered (ferromagnetic) state with {M) 7^ 0. The different phases of the material are thus 
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conveniently characterized by the average magnetization, which for this reason is known as 
the order parameter of the ferromagnetic phase transition. 

In the phenomenological approach put forward by Landau the above mentioned temper- 
ature dependence of the equilibrium order parameter (M) is reproduced by anticipating 
that the free-energy density of the system at a fixed but not necessarily equilibrium value 
of the average magnetization has the following expansion 

f((M),T) ~ /(0,T) + a(T)(M) 2 + ^-(M)' (4) 

for small values of (M), and that the coefficients of this expansion behave near the critical 
temperature as 

a(T) ^a (^- l) (5) 



and 



(3{T) ~ p , (6) 

respectively, with a and j3 positive constants. 

Hence, above the critical temperature a(T) and f3{T) are both positive. As a result the 
free energy is minimal for (M) = 0, which corresponds exactly to the paramagnetic phase. 
Moreover, for temperatures below the critical one at(T) is negative and the free energy is 



indeed minimized by a nonzero average magnetization with magnitude J—a(T)/j3(T). Just 
below the critical temperature the latter equals 



'*» V A V T, 

which after substitution in Eq. (f|) gives rise to an equilibrium free-energy density of 



f((M) eg ,T)^f(0,T)-^-(l--) . (8) 



Therefore, the second derivative d 2 f/dT 2 is discontinuous at the critical temperature and 
the phase transition is of second order according to the Ehrenfest nomenclature. 



Note that minimizing the free energy only fixes the magnitude and not the direction of 
(M) . This degeneracy is caused by the fact that the Hamiltonian in Eq. (Q) is symmetric 
under an arbitrary rotation of all the spins Si. Consequently, the free energy must be 
symmetric under a rotation of the average magnetization and only even powers of (M) can 
appear in its expansion (cf. Eq. (^)). Due to this behavior the ferromagnet is a good example 
of a system with a spontaneously broken symmetry, i.e. although the Hamiltonian is invariant 
under the operations of a group, its ground state is not. In the case of a ferromagnet the 
symmetry group is 5*0(3), which is broken spontaneously below the critical temperature 
because the average magnetization points in a certain direction. Which direction is chosen 
in practice, depends on the surroundings of the system and in particular on (arbitrary small) 
external magnetic fields that favor a specific direction. 

After this summary of the Landau theory we are now in a position to introduce two 
time scales which turn out to be of great importance for the nucleation of Bose-Einstein 
condensation. To do so we consider the following experiment: Imagine that we have a piece 
of ferromagnetic material at some temperature T\ above the critical temperature. Being in 
thermal equilibrium the material is in the paramagnetic phase with (M) = 0. We then 
quickly cool the material to a new temperature T 2 below the critical temperature. If this is 
done sufficiently fast, the spins will have no time to react and we obtain a nonequilibrium 
situation in which the free energy has developed a 'double-well' structure but the average 
magnetization is still zero. This is depicted in Fig. |](a). In such a situation there is a 
typical time scale for the relaxation of the average magnetization to its new equilibrium 
value \J— a(T 2 )//?(T 2 ), which we denote r co h- 

However, in the case of magnetically trapped atomic hydrogen, the gas is isolated from 
its surroundings and it is not possible to perform the cooling stage mentioned above. As a 
result the gas has to develop the instability associated with the phase transition by itself. 
The time scale corresponding to this process is called T nuc i and is schematically shown in 
Fig. ^](b). Combining the two processes we are lead to the following physical picture for the 
nucleation of Bose-Einstein condensation. After the quench into the critical region the gas 
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develops an instability on the time scale r nuc i. On this time scale the actual nucleation takes 
place and a small nucleus of the condensate is formed, which then grows on the time scale 
r co h as shown in Fig. |2|. To solve the nucleation problem we are thus left with the actual 
determination of these two time scales. Clearly, before this can be done we need to know 
the correct order parameter of the phase transition. 



B. Order parameter 



Ever since the pioneering work of Bogoliubov |L8[ it is well known that the order param- 
eter for Bose-Einstein condensation in a weakly-interacting Bose gas is a somewhat abstract 
quantity, which is most conveniently discussed by using the method of second quantization. 
In this method all many-body observables are expressed in terms of the creation and annihi- 
lation operators of a particle at position x denoted by if)^(x) and ip(x), respectively |0|. For 
example, for a gas of particles with mass m and a two-body interaction potential V(x — x') 
the Hamiltonian equals 

H = f dxip\x) ~ h V ip{x) + - fdx [ i) ] {x)i)\x')V{x- x')i){x')i){x) (9) 



2m v ' 2 
and the total number of particles is given by 



N — dx^{x)ip(x) . (10) 



The method is also particularly useful for a Bose system because the permutation symmetry 
of the many-body wave function is automatically accounted for by assuming the commuta- 
tion relations [i/i(x), ij)(x*)\ = [ip^(x), ip\x')} = and [ip(x), ip\x')} = 8(x — x') between the 
creation and annihilation operators. 

In the language of second quantization the order parameter for the dilute Bose gas is 
the expectation value (ip(x)). Analogous to the case of the ferromagnetic phase transition, 
a nonzero value of this order parameter signals a spontaneously broken symmetry. Here the 
appropriate symmetry group is U(l), since the Hamiltonian of Eq. (|Sp is invariant under the 
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transformation ip{x) — > xjj(x)e 1 ^ and ij)^(x) — > ^ (x)e _M? of the field operators, whereas their 
expectation values are clearly not. Notice that the U(l) symmetry of the Bose gas is closely 
related to the conservation of particle number. This is most easily seen by observing that 
the invariance of the Hamiltonian is due to the fact that each term in the right-hand side of 
Eq. (|9|) contains an equal number of creation and annihilation operators. The relationship 
can also be established in a more formal way by noting that the U(l) gauge transformations 
are generated by the particle number operator. As we will see later on, it has important 
consequences for the dynamics of the order parameter. 

To understand why (ip{x)) is the order parameter associated with Bose-Einstein conden- 
sation, it is convenient to use a momentum-space description and to introduce the annihila- 
tion operator for a particle with momentum hk 

H = \ dx^{x)—^=- (11) 

and the corresponding creation operator at' by Hermitian conjugation. The basis of states 
for the gas is then characterized by the occupation numbers {-/VV}. If the gas is condensed, 
there is a macroscopic occupation of the zero-momentum state and the relevant states are 
\Ntf, {Nfi,}^^) with only Nq proportional to N. Within this subspace of states we have 

(offtflj) = (N S ) * (N 5 ) + 1 = (ogofft) (12) 
and we can neglect that and do not commute. As a result we can treat these operators 



as complex numbers |18] and say that (a^a^) = (oq )(°o) or equivalently that (aj) = yiVg. 
In coordinate space the latter reads (i/j(x)) = y^j, with = N^/V the condensate density. 

The above argument essentially tells us that a sufficient condition for a nonzero value of 
(ip(x)) is (Nq) ^> 1. Although this is intuitively appealing, it is important to point out that 



it is not generally true. Consider for example the ideal Bose gas [ 2U| . In the grand canonical 
ensemble the total number of particles in the gas is given by 

k k s 



where (3 is l/k B T, eg is the kinetic energy h 2 k 2 /2m, £ is the fugacity and /i is the 
chemical potential. 

At high temperatures the fugacity is small and we are allowed to take the continuum 
limit of Eq. (fl3|) , which results in the equation of state 

n = ^3/ 2 (C) , (14) 



using the thermal de Broglie wavelength A = J 2tcTi 2 /mk B T and the Bose functions g n (() 
defined by 

^(0 = ^/ cfa f, ■ (15) 
1 [n) Jo C e — 1 

Lowering the temperature while keeping the density fixed, the fugacity increases until it 
ultimately reaches the value one at the critical temperature 

_ 2nh 2 f n \ 2/ \ 2nh 2 f n \ 2 / 3 
°~ mk B {g 3/2 (l)J ~ mk B \2M2) ' ( } 

At this point Eq. (|14|) ceases to be valid because the occupation number of the zero- 
momentum state, which is equal to — £), diverges and must be taken out of the discrete 
sum in Eq. (|13D before we take the continuum limit. Moreover, we only need to treat the 
zero-momentum term separately because in the thermodynamic limit the chemical poten- 
tial goes to zero as V^ 1 , whereas the kinetic energy for the smallest nonzero momentum 
decreases only as V~ 2 ^ 3 . Consequently, below the critical temperature the equation of state 
becomes 

n = n o + ^3 03/2(1) (17) 

and leads to a condensate density equal to 

"^"t 1 "© 3 ' 2 )- (18) 

We thus find that the average occupation number (N$) is at all temperatures given by 
£/ (1 — £), i- e - hs value in the grand canonical ensemble with the density matrix e~^ H-MAr ). 
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Since this density matrix commutes with the particle number operator, we conclude that 
in the case of an ideal Bose gas there is a macroscopic occupation of the zero-momentum 
state without a spontaneous breaking of the U(l) symmetry. To show more rigorously that 
{ip{x)) = at all temperatures we determine the free-energy density of the gas as a function 
of the order parameter (ip(x)). Dealing with a noninteracting system it is not difficult to 
obtain 

m(x)),T) = -^T)mx)}\ 2 (19) 

for a homogeneous value of the order parameter. Because fi < the minimum is indeed 
always at (4>(x)) = and it is necessary to identify the condensate density with the order 
parameter of the ideal Bose gas (cf. Eq. (pT8[)). 

Notwithstanding the previous remarks, the order parameter for Bose-Einstein condensa- 
tion in a weakly-interacting Bose gas is given by (ip(x)). This was put on a firm theoretical 



basis by Hugenholtz and Pines ||21|| , who calculated the free energy as a function of the 
above order parameter and showed that at sufficiently low temperatures the system devel- 
ops an instability that is removed by a nonzero value of {ip(x)). In addition, they derived an 
exact relationship between the chemical potential and the condensate density, which turns 
out to be valid also in the nonequilibrium problem of interest here and is important for an 
understanding of how the U(l) symmetry is broken dynamically. 



C. Condensation time 

We have argued that by means of evaporative cooling a doubly-polarized atomic hydrogen 
gas can be quenched into the critical region of the phase transition and that this kinetic 
part of the condensation process is described by a quantum Boltzmann equation. As a result 
the gas acquires on the time scale r e \ an equilibrium distribution with some temperature T, 
which is slightly above the critical temperature To of the ideal Bose gas because a condensate 
cannot be formed at this stage. 



11 



For the study of the subsequent coherent part of the condensation process it is therefore 
physically reasonable to assume that at a time t the density matrix p(to) of the gas is well 
approximated by the density matrix of an ideal Bose gas with temperature T. The evolution 
of the order parameter (ip{x)) for times larger than to is then completely determined by the 
Heisenberg equation of motion 

ih^^ = mt),H], (20) 

for the field operator. Substituting herein the Hamiltonian of Eq. (§) and taking the expec- 
tation value with respect to p(t ) we find 

ih m 2 t)} = ^^(^(x,t)) + J dx*V{x-x'){^\x\t)^,t)^{x,t)) , (21) 

where the complicated part is of course the evaluation of (ift (a?, t)ip(x', t)ijj(x, t)). In lowest 
order we simply have 

(ij\x',t)ij(x',t)ij(x,t)) ~ (^(f ,t))(il>(x",t)) (i/;(x,t)) + (^(x',t)^(x',t))(i/;(x,t)) 

+ (^(f / ,t)^(f,t))^(f,t)), (22) 

which after substitution into Eq. (plf) leads to 

ih jt + ^r) ' t))= I dS ' v{s ~ s ' ] ( ^ t(f '' tMf/ ' t]) ' t]) 

+ J dx'V(x-x')(^(x',t)%l)(x, t)}(iP(x',t)} 
+ J dx* V(x-x')(^(x',t))(^(x',t))(^(x,t)) , (23) 

and thus corresponds exactly to the Hartree-Fock approximation. 

To proceed we must restrict ourselves to the case of a dilute Bose gas in the quantum 
regime. Introducing the scattering length a, which is of the order of the range of the 
interaction, the quantum regime is characterized by a/A 1. We therefore need to consider 
only s-wave scattering and can neglect the momentum dependence of various collisional 
quantities. In particular, we can replace the potential V(x — x') by the contact interaction 
V^5(x — x') with Vq = J dxV(x). Hence, in the Hartree-Fock approximation we obtain 

12 



{ ih Jt + { ^ t]) = ( 2nV o+Vo\mt,t))\ 2 ) (1>frt)) , (24) 

having only the trivial solution (i[)(x, t)) = for a space and time-independent order param- 
eter. Within this lowest order approximation we thus conclude that T nuc \ = oo and that the 
formation of a condensate will not take place. 

Fortunately, it is well known that the Hartree-Fock approximation is not sufficiently ac- 
curate for a dilute Bose gas because the diluteness condition na 3 <C 1 implies that we should 
consider all two-body processes, i.e. two particles must be allowed to interact also more than 
once. The appropriate approximation is therefore the ladder or T-matrix approximation and 
is diagrammatically explained in Fig. |^. Moreover, in the degenerate regime where the tem- 
perature T is slightly larger than To and the degeneracy parameter nA 3 is of 0(1), the 
condition a/A ^ 1 implies that also naA 2 <C 1 or physically that the average kinetic energy 
of the gas is much larger than the typical interaction energy. Consequently, an accurate 
discussion of the nucleation of Bose-Einstein condensation in a weakly-interacting Bose gas 
requires an evaluation of (^(x 1 , t)ip(x', t)ip(x, t)) within the T-matrix approximation and in 
zeroth order in the gas parameters a/ A and naA 2 . 

Although it is easy to formulate this objective, to actually perform the calculation is 
considerably more difficult. It is most conveniently accomplished by making use of the 



Keldysh formalism |22] which has been reviewed by Danielewicz |23| using operator methods. 
For a functional formulation of this nonequilibrium theory and for the technical details of the 
somewhat tedious mathematics we refer to our previous papers |[L6|| . Here we only present 
the final results and concentrate on the physics involved. 

Due to the fact that we are allowed to neglect the (relative) momentum dependence of 
the T matrix, the equation of motion for the order parameter (ijj(x,t)) acquires the local 
form of a time-dependent Landau- Ginzburg theory 

{ ih jt + ^r) Wf,t)) = {s {+) (t)+T^\mx,t)}\ 2 ) (v>(*,*)> , (25) 

which is recovered from a variational principle if we use the action 

13 



r r ( d f> 2 X7 2 T(+) \ 

S((mt)),T) = jdtj dx^&t))* [ifi- + - fif(+)(t) - — |(^(x,t))| 2 J tyfrt)) . 

(26) 

Here S^ + \t)5(t — t') is a good approximation for the retarded self-energy hT,^ + \0; t, t') of 
a hydrogen atom with zero momentum and T"W ~ Aith 2 a/m is the effective interaction 
between two such atoms. Clearly, the action in Eq. (|26|) is the desired generalization of the 
Landau free energy and corresponds precisely to the physical picture presented previously 
in Figs. |(b) and|(b). 

Therefore, r nuc i is determined by the time dependence of the coefficient which is 
shown in Fig. f| for three different initial temperatures. If the temperature T is much 
larger than T , S^ + '(t) is constant and equal to Sufi 2 an /m. In this region of the phase 
diagram coherent processes are negligible and the evolution of the gas is described by a 
Boltzmann equation. Lowering the temperature, the occupation numbers for momenta 
hk < 0(h/A) rise and lead to an enhancement of the coherent population of states with 
momenta hk < 0(hy/na) 0(h/A). This is signaled by the increasing correlation length 
£ = fij y / 2mS'( + )(oo). At the critical temperature T c = T (l + O(a/A )) we have S^ + \oo) = 
and the correlation length diverges. Below that temperature, but still above To so as not 
to have a condensate already in the initial state, S^ + '(t) actually changes sign and the gas 
develops the required instability for a Bose-Einstein condensation. The change of sign takes 
place at 

(s( - ( » + (tM7^T))' (27) 

which shows that r nuc i is in general of 0{fi/kBT c ) except for temperatures very close to the 
critical temperature. Clearly, this time scale is due to the fact that all states with momenta 
hk < 0(h/A) cooperate in the coherent population of the one-particle ground state. 

After a small nucleus of 0(n(a/A c ) 2 ) has been formed, the subsequent buildup of the 
condensate density is determined by the equation of motion Eq. (P5|). Looking at the right- 
hand side we immediately see that the time scale r co h involved in this process is typically 
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of 0{h/n^T^ + >) or equivalently of 0((h / kBT c )(l / n$ak 2 c )) . Therefore, r co h ^> r nuc i as antic- 
ipated in Fig. 0. The physical reason for this time scale is that after the nucleation of the 
phase transition the buildup of the condensate density is accompanied by a depopulation of 
the momentum states with hk < 0(Ti^Jn^a). As a result it is not difficult to show that in the 
limit t — > oo the condensate density is of 0(na/A c ) and thus that T co h = 0((A c /a) 2 h/kBT c ). 

Finally, it is interesting to point out how the gas can conserve the total number of 
particles and apparently at the same time break the U(l) gauge symmetry that is responsible 
for this conservation law. To that end we write the field operator ip(x,t) as the sum of its 
expectation value {ip(x,t)) and the fluctuation tp'(x,t), and introduce a time-dependent 
chemical potential fi(t) by means of 

(if)(£,t)) = sf^exp^-^^dt'^t')^ . (28) 



Substituting the latter in the action of Eq. (|2q) and minimizing with respect to yn^t) gives 
for t > t c 

which determines the growth of the condensate density and is in effect a nonequilibrium 
version of the Hugenholtz-Pines theorem ||21||. Furthermore, by considering the fluctuations 



around \ln^{t) we can show that the chemical potential is determined by the constraint 

n = n 5 (t) + ^ Jdx (ip'^x, t)i/>'(x, t)) , (30) 

enforcing the conservation of particle number at all times. In the complex plane (ijj(x,t)) 
thus moves radially outward along a spiral as shown in Fig. [5]. Consequently, the phase of 
the order parameter has never a fixed value and the U(l) symmetry is not really broken 
dynamically. This is of course expected since the system evolves according to a symmetric 
Hamiltonian. 



15 



III. CONCLUSIONS AND DISCUSSION 



We studied the evolution of a doubly-polarized atomic hydrogen gas in a magnetic trap 
and showed that by means of evaporative cooling the gas can accomplish the Bose-Einstein 
phase transition within its lifetime T; ne ;. The condensation process proceeds under these 
conditions in three stages: In the first kinetic stage the gas is quenched into the critical region 
Tq < T < T c . A typical time scale in this part of the evolution is given by the time between 
elastic collisions r e i, which for a degenerate gas is of 0((A c /a) 2 h/kBT c ). In the following 
coherent stage the actual nucleation takes place on the time scale r nud = 0{h/k,BT c ) by 
means of a coherent population of the zero-momentum state. The small nucleus formed 
in this manner then grows on the much longer time scale T co h = 0((A c /a) 2 h/kBT c ) by 
a depopulation of the low-momentum states, having hk < 0(h^n^a). In the third and 
last stage of the evolution the Bogoliubov quasiparticles produced in the previous stage 
have to come into equilibrium with the condensate. This process can again be treated 



by a kinetic equation and was studied by Eckern p4|, who found that the corresponding 
relaxation time r re \ is of 0((A c /a) 3 h/kBT c ). In the case of atomic hydrogen this turns out 
to be comparable to the lifetime of the system. Summarizing, we thus have the sequence 
i~nud T C oh — T ei ^ r rei — T inei for the various time scales involved in the phase transition. 
The most important requirement for the achievement of the phase transition is therefore 
T ei ^ T ineh which is relatively mild and should not pose an insurmountable problem for 
future experiments aimed at the realization of Bose-Einstein condensation. 

Having arrived at this conclusion, it is necessary to discuss a recent paper by Kagan, 



Svistunov, and Shlyapnikov [25| that also considers the evolution of a weakly-interacting 
Bose gas after the removal of the most energetic atoms. In this paper the authors agree that 
the evolution of the gas is divided into a kinetic and a subsequent coherent stage. Moreover, 
their detailed study of the kinetic part of the evolution confirms our conjecture that the gas 
is quenched into the critical region on the time scale T e \. The investigation of the coherent 
part, however, leads to the extreme result that a Bose-Einstein condensation cannot occur 
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in a finite amount of time. To understand why this conclusion is reached we briefly present 
their line of thought. 

At the end of the kinetic stage the gas has acquired large average occupation numbers 
for the states with momenta hk < hko = O(hy/noa), where no is the density of particles 
with these small momenta. Therefore, Kagan, Svistunov, and Shlyapnikov argue that for a 
study of the coherent part of the evolution we must use the initial condition 

„ik-x 



(#Mo)>= E (31) 

k<k V v 

together with the nonlinear Schrodinger equation 

in m^ L = ^^ +Ti+)mm ^ {mt)h (32) 

which has the equilibrium solution (ip(x,t)) = y/noexp(—ifx t) and /z = n T^ + \ Conse- 
quently, all the particles that have initially momenta hk < hko are in the limit t — > oo 
assumed to be in the condensate. 

Linearizing the Hamiltonian of Eq. Q3"2]) around this equilibrium solution they then 
observe that the energy involved with a magnitude fluctuation of the order parameter is 
+ n T^ + \ whereas the energy involved with a phase fluctuation is only eg. As a result 
they assert that on the time scale T amp i = r co h = 0(h/noT^) a state is formed in which the 
amplitude of (ip(x,t)) is fixed, but the phase is still strongly fluctuating because the cor- 
responding time scale r p h is much longer and even diverges as V 2 ^ 3 in the thermodynamic 
limit. Hence, for finite times the gas is in a state with a so-called quasicondensate |26[ and 
a real condensate is only formed in the limit t — > oo. 

Clearly, this physical picture of two different time scales for the amplitude and phase 
fluctuations of the order parameter is only applicable if these fluctuations exist indepen- 
dently of each other. Looking only at the Hamiltonian this indeed seems to be the case. 
However, a correct discussion of the fluctuations must be based on the equations of mo- 
tion or equivalently the Lagrangian. The latter contains a first-order time derivative which 
strongly couples the amplitude and phase fluctuations. Therefore, a dilute Bose gas does 
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not have two but only one dispersion relation, i.e. the well-known Bogoliubov dispersion 



e fe( e fe + 2hoT"W), and we are lead to r ph = r amp i. It is interesting to note that in the case 
of a neutral BCS-type superfluid we do have two different time scales because the Lagrangian 
now contains a second-order time derivative and the amplitude and phase fluctuations are 
indeed independent in lowest order PT| , pE| . 

An even more serious problem with the approach of Kagan, Svistunov, and Shlyapnikov 
is their claim that the use of the initial condition in Eq. (|3lD is justified because (N^) 3> 1. 
As we have pointed out before this is not true in general. For (ip(x, t)) to be nonzero we must 
show that the system has a corresponding instability. However, within the T-matrix approx- 
imation it is not difficult to show that the instability associated with a quasicondensate is 
always preceded by the instability corresponding to the formation of a condensate. This 
implies that we always have to take Bose-Einstein condensation into account first. After 
that has been accomplished by means of the theory reviewed here, it is of course no longer 
relevant to consider the appearance of a quasicondensate. 
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FIGURES 

FIG. 1. Visualization of (a) the time scale r co h for the relaxation of the order parameter to its 
equilibrium value and (b) the time scale T nuc i associated with the appearance of the instability. 

FIG. 2. Visualization of the time scales r co h and T nuc i, using the time dependence of (a) the 
order parameter and (b) the coefficient a of the quadratic term in the free energy. 

FIG. 3. Diagrammatic representation of the T-matrix equation. The wavy line corresponds to 
the interaction and the straight line to the non-interacting one-particle Green's function. 

FIG. 4. Time dependence of the coefficient for three different initial temperatures of the 
Bose gas. 

FIG. 5. Evolution of the complex order parameter (tjj(x)), which is constrained by the require- 
ment of particle number conservation. 
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